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ABSTRACT 

Galaxy-scale strong gravitational lenses with measured stellar velocity dispersions allow a test of the weak- 
field metric on kiloparsec scales and a geometric measurement of the cosmological distance-redshift relation, 
provided that the mass-dynamical structure of the lensing galaxies can be independently constrained to a suf- 
ficient degree. We combine data on 53 galaxy-scale strong lenses from the Sloan Lens ACS Survey with a 
well-motivated fiducial set of lens-galaxy parameters to find (1) a constraint on the post-Newtonian parameter 
7 = 1.01 ±0.05, and (2) a determination of fl\ = 0.75 ±0.17 under the assumption of a flat universe. These 
constraints assume that the underlying observations and priors are free of systematic error. We evaluate the 
sensitivity of these results to systematic uncertainties in (1) total mass-profile shape, (2) velocity anisotropy, 
(3) light-profile shape, and (4) stellar velocity dispersion. Based on these sensitivities, we conclude that while 
such strong-lens samples can in principle provide an important tool for testing general relativity and cosmology, 
they are unlikely to yield precision measurements of 7 and fl\ unless the properties of the lensing galaxies are 
independently constrained with substantially greater accuracy than at present. 

Subject headings: cosmological parameters — gravitational lensing — relativity 



1. INTRODUCTION 

Einstein's theory of general relativity (GR) has been an ex- 
tremely successful description of gravity. It has passed all 
current experimental tests, most famously Eddington's mea- 
surement of light d eflection during the solar eclipse of 1919 
jDyson et al.lll92.0l); observ ation of the gravitational redshift 
by iPound & Reb ka (I960]); the successfu l operation of the 
Global Positionin g Satellites (IAs ibvl2002l); measure ments of 
the Shapiro delay dShapirol 1964 iBertotti et alJ2003l) : and ex- 
tensive studies of relativistic effects in binary radio pulsar sys- 
tems, including verification of energy loss via gravitational 
wave s as observed in the Hulse-Taylor pulsar dTavlor et al.l 
119791) . Tests of gravity at ever higher precisions continue 
to be pursued through techniques such as lunar laser rang- 
ing, where the Earth- Moon separation is pr ecisely measured 
as a function of time (Williams et al. 2 0041). The parameter- 
ized post-Newtonian (PPN) framework dThorne & Willi 1 97 lh 
provides a systematic, quantitative way in which to formulate 
and interpret tests of gravity. The post-Newtonian parame- 
ter traditionally denoted by 7 — a measure of the amount of 
spatial curvature per unit mass — is currently c onstrained to be 
7= 1 + (2.1 ±2.3) x 10" 5 (iBertotti et al.l2003l) on solar-system 
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length scales. Within the post-Newtonian parameterization of 
GR, 7 has no scale dependence, and thus constraints on de- 
viations from 7 = 1 on significantly larger length scales can 
provide further tests of the theory. 

GR has had further success in its application to physical 
cosmology. From assumptions of homogeneity and isotropy, 
the evolution of the universe as a function of time can be pre- 
dicted from the knowledge of the densities of its constituents. 
While it was initially assumed that the density of the universe 
was dominated by matter (including a significant component 
of non-baryonic "dark matter"), the discovery that the uni- 
verse is expanding at an accelerating rate (as trace d by type la 
supernovae) has forced a mod ification of this idea (Ri ess et all 
[1991 iPerlmutter et al.|[l999l) . It is now thought that the uni- 
verse is filled with a "dark energy", an unknown substance 
with negative pressure. The composition of the universe de- 
termines its evolutionary history through the FRW metric and 
the Friedmann equations, and the evolution history of the uni- 
verse in turn sets the distance-redshift relation. Thus, mea- 
surements of the conversion from redshift to distance con- 
strain the density of dark energy, which we denote simply by 
SI A in this work (i.e., assuming an equation of state P = -pc 2 ). 

This paper presents an analysis of the quantitative con- 
straints that can be placed on GR as a theory of gravity 
and on the density of dark energy in the universe using the 
large and homogeneous set of strong gravitational lens galax- 
ies observe d by the Sloan Len s ACS (SLACS) Survey col- 
laboration (Bol ton et alj|2006al) . To date, research using the 
SLACS sample has been primarily focused on constraining 
the internal structure and dynamics of the SLACS lenses un- 
der the assump t ion of GR and a ACDM cosmology (e.g . 
Treu et alj|2006t iKoopmans et al.|[2006t iGavazzi et al.fl2007l: 
Koop mans et al.l 120091) . However, with reasonable prior as- 
sumptions and independent measurements of the structure 
of early-type galaxies such as the SLACS lenses, the prob- 
lem can be inverted and the sample used to constrain the 
parameters of GR and of the cosmology. The first avenue 
was explored for an initial sample of 15 SLACS lenses by 
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iBolton et al.l d2006bl) . who found the post-Newtonian param- 
eter 7 to be 0.98 ± 0.07 on kiloparsec scales by adopting 
priors on early-type galaxy structure taken from observa- 
tions of the local univer se. The second av enue, originally 
suggested in passing by iGolse et al.l (|2002), was examined 
at length by iGrillo et al.l (120081) . who discussed the use of 
the SLACS lensing systems as a means to determine cos- 
mological parameters, finding that the "concordance" val- 
ues of (^ m ,^a) ^ (0.3,0.7) fell within their 99% confidence 
regions. In both cases, the general approach consists of 
using angular "Einstein radii" measured from imaging, in 
combination with spectroscopic redshifts and parameterized 
galaxy mass models, to predict an observable stellar velocity 
dispersion — which explicitly depends upon the parameter of 
interest, whether post-Newtonian or cosmological — and then 
comparing this predicted value with the velocity dispersion 
measured from spectroscopy. Alternatively, these techniques 
can be characterized as treating strong-lensing Einstein radii 
as "standard angles" calibrated by the s tellar velocity disper- 
sions of the lensing galaxies. (Also see Gavazzi et al.l (120081) 
for a related cosmological application.) 

In this work we re-examine the constraints that can be set 
on both 7 and Qj\ using strong gravitational lensing by early- 
type galaxies, given the availabili t y of the expanded SLACS 
sample presented by Bolton et al.l d2008al) . Our final empha- 
sis will not be as much on the most probable values for these 
parameters, but rather on the uncertainties due to systematic 
errors made in estimating the properties of the lensing galax- 
ies. In Section 12 we discuss the SLACS lens sample that is 
used in our analysis. The lensing formalism that we utilize 
is reviewed in Section [3] Our determination of 7 and SI a for 
a fiducial set of lensing parameters is presented in Section [4] 
along with a brief discussion of the statistical analysis that was 
employed. Section is devoted to an in-depth evaluation of 
the effects that systematic errors in the lens parameters have 
on the nominal values for 7 and SIa- Finally, in Section|6]we 
summarize our findings and draw some conclusions on how 
to improve the accuracy of strong lensing results in the future. 

2. SLACS LENS SAMPLE 

The post-Newtonian strong-lensing constraints published 
by IBolton et"al] d2006bl) were based on a sample of 15 strong 
galaxy-galaxy gravitational lenses from the S L ACS Sur- 
vey, w hich are described in IBolton et ail d2006al) . iTreu et ail 
d2006l) . and lKoopmans et al.l (12006b . The current work uses a 
more recent and larger sample of strong lenses from this sam e 
survey, which is described in detail by IBolton et al.l (I2008al) . 
These systems were all selected from within the sp ectroscopic 
datab ase of the Sloan Digital Sky Survey (SDSS: lYork et al.l 
120001) based upon the presence of two significantly different 
galaxy redshifts within a single spectrum, obtained with a 3"- 
diameter fiber aperture. As a consequence of explicit lensing 
and other selection effects, the SLACS sample consists pri- 
marily of massive, early-type (i.e., elliptical and SO) galaxies 
lensing much fainter and more distant emission-line galaxies. 

For the present work, the key observables in each system 
are the redshifts of the two components (foreground "lens" 
and background "source"), the stellar velocity dispersion of 
the lens galaxy, a power-law index representing the charac- 
teristic slope of the luminosity profile, and the angular Ein- 
stein radius of the strongly lensed image of the more distant 
galaxy. The first three of these quantities are measured from 
SDSS spectroscopy, while the last two are measured from 
high-resolution follow-up images obtained by the SLACS sur- 



vey through the F814W (/-band) filter with the Wide Field 
Channel of the Advanced Camera for Surveys (ACS) aboard 
the Hubble Space Telescope (HST). Note that mass and light 
models fitted to the HST data include a projected axis-ratio 
parameter, q, to capture ellipticity. To connect these mod- 
els to the axisymmetric approximation of the analytic Jeans 
equation-based framework introduced below, we use the in- 
terchange 



R 



R q = \l qx 2 +y 2 /q 



(1) 



which conserves the total mass or light within a given iso- 
density or isobrightness contour and is consistent with the 
Einstei n-radius and effectiv e-radius measurement conven- 
tions of IBolton etaD (l2008al) . 

The 53 lenses co nsidered in this work are those from 
IBolton et al.ld200 8a) that have single lens galaxies with early- 
type morphology, successful quantitative strong lens models 
(and hence measured Einstein radii), and sufficiently high 
signal-to-noise measurements of lens-galaxy stellar veloc- 
ity dispersion. These lens galaxies have redshifts in the 
range z ~0. 1-0.3, and given these relatively low redshifts 
(and hence modest look-back times), there is unlikely to 
be significant structural difference between the SLACS lens 
galaxies and the loca l-universe elliptical galaxies that we and 
IBolton et al.l (l2006bl) use as a calibration sample. 

3. LENSING FORMALISM 

We start with the weak-field PPN form of the Schwarzschild 
metric 



dr 2 =dt 2 {\-2M/r)-dr 2 (\-2 1 M/r)-r 



(2) 



with G = c= 1 , where 7 is a parameter to be constrained, and 
7=1 for GR. For thin gravitational lenses, the "lens equation" 
which governs the deflection of light is then 



(l+ 7 ) 



vv>(0) , 



(3) 



where 6 S is the angular location of the source, and 9 is the an- 
gular location of the image. The scaled, projected Newtonian 
potential tp is defined as 



D 



LS 



D L D S c 2 



dZ $(D L 8,Z) 



(4) 



where $ is the traditional Newtonian potential and Z is the 
line-of-sight distance between the observer and the source, 
with Z = at the lens. Here, D s is the distance to the source, 
Dl is the distance to the lens, and Dig is the distance between 
the lens and the source; all three are cosmological angular- 
diameter distances. 
For a lens that is cylindrically symmetric around the line 

of sight, the Einstein ring radius is found by setting 6 S = to 
yield 

(1+7) (dm\ 
2 I de ■ 



Or 



(5) 



The Einstein radius provides the characteristic angle (or 
length) in gravitational lensing. The angular size of the Ein- 
stein radius corresponding to a point mass M is given by 



-m 



4GM D LS 

C 2 DsDr 



1/2 



(6) 
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For more general cylindrically symmetric mass distributions 
with respect to the line-of-sight, this continues to hold with 
M = Me = M(Re), the mass contained within a cylinder of 
radius equal to the Einstein radius. Rearranging terms, and 
noting that R E = D L 9 E , we find: 



c 2 D s 



GM E 



Re (1+7)4D ls 



(7) 



a form which will prove useful. Regardless of the extent of the 
mass distribution, only the mass interior to the Einstein radius 
has a net effect on the deflection of the ligh t in the circularl y 
symmetric (or, more generally, homoeoidal dSchrammll 1 9901) ) 
case. 

We will mainly consider a flat cosmology, where £1^ = 
and thus Om + Oa = 1 ■ This means there is only one free cos- 
mological parameter appearing in the distance-redshift rela- 
tion. This restriction on the dimension of the parameter space 
yields a sharper constraint (see Section [?}. The assumption 
of flatness is not an unreasonable one. Curr ently, indepen- 
dent m easurements show |1— Om _ Oa| < 02 (iKomatsu et al.l 
120091). and t he theoretical prediction of inflation is that O ss 1 
dGuMl98l . 

In this flat case, the relation among the angular diameter 
distance between two objects {Dam) and their line-of-sight 
comoving distances Dq is 



D A „ = 



1 



1+Z2 



(D c2 -D cl ) 



(8) 



where z is the redshift. Then, the ratio of angular diameter dis 
tances Dls /Ds which appears in our expression for the Ein 
stein radius is simply 

Dls_ 1 Dcl 
D s D cs 



(9) 



Since this is a ratio of two comoving distances, there is no 
dependence on the Hubble constant. Thus, this ratio depends 
only on the source and lens redshifts, which are measured, 
and the choice of OA- 
One of the simplest models one can write down for an ellip- 
tical galaxy is a scale-free model based on power-law density 
profiles for the total mass density, p, and luminosity density, 



(10) 



p{r) = Po[ — 
Jo 



v{r) = Vq\ — 
Jo 



(11) 



(see, e.g.. lKoopmansl2006l) . In general, the three-dimensional 
velocity dispersion tensor is not isotropic. The deviation can 
be characterized through the anisotropy parameter (3, defined 

as 

P(r) = l-af/af. (12) 

where a r is the radial value. Having assumed spherical sym- 
?, and we refer to this component as "tan- 



metry, a? = a\ 



gential". For the purpose of the current analysis it will suffice 
to assume that is indep endent of r. This is the model used b y 
by iBolton et all (l2006bl) and discussed by iKoopmansI (120061) . 

We define r to be the spherical radial coordinate from the 
lens center and R to be the cylindrical radius, i.e., perpendicu- 
lar to the line of sight (defined to be along the Z-axis). So by 
definition r 2 = R 2 + Z 2 . 



A key ingredient in all the lensing measurements is the ob- 
served velocity dispersion, which is a projected, luminosity 
weighted average of the radially-dependent velocity disper- 
sion profile of the lensing galaxy. In order to predict this value 
based on a set of galaxy parameters, we start with an expres- 
sion, derived from the spherical Jeans equation, for a 2 (r), the 
radial velocity dispersion of the luminous matter (i.e., stars) 
(lBinnevHl980f) : 



. G dr' v{r')M{r'){r') 2 P- 
a r {r)= '■ 



V{f) 



(13) 



for the case where the velocity anisotropy parameter, j3, is a 
constant. Note that the use of Equation $13[ is based on the 
assumption that the relationship between stellar number den- 
sity and stellar luminosity density is spatially constant, an as- 
sumption unlikely to be violated appreciably within the effec- 
tive radius of the early-type lens galaxies under consideration. 
Using the mass density profile in Equation ( [Tol l, it is straight- 
forward to show that the relation between the mass contained 
within a spherical radius r and M E is 



M(r): 



0rA(a) \R E 



3-Q 



M E 



(14) 



where we have defined the ratio of gamma functions 

A(x) = r(^)/r(|) . 

Therefore, after evaluating the integral in Equation ( fT3l we 
find 



(15) 



GM E 



Re 



Rr 



(16) 



0F(£-2/3)A(a) 

where, following IKoopmansI ( 120061) . we have defined £ = 
S + a-2 as a convenient combination of the power-law ex- 
ponents. 

Because all of the elliptical galaxy lenses in this study 
have been directly imaged with HST, we can measure their 
projected two-dimen sional luminosity profiles. Following 
IBolton et al.l d200 6b). we directly fit point-spread function 
(PSF) convolved two-dimensional power-law ellipsoid im- 
ages to HST F814W imaging data over a circle of radius l."8 
centered on the lens galaxies, adding 1 to the best-fit two- 
dimensional power-law index so as to account for a deprojec- 
tion into three dimensions to obtain our adopted 5 values for 
each system. 

The velocity dispersions are measured from SDSS spectra. 
Each SDSS spectroscopic fiber subtends a circle of radius 1 "5 
on the sky, with atmospheric blurring adding an additional 
consideration. The actual aperture weighting function w(R) 
should thus be the convolution of the atmospheric seeing, rep- 
resented by a Gaussian, 



s(R) = exp 
and the fiber aperture 

a{R) = n 



R 2 
2cr 2 



R 



(17) 



(18) 



where n is the rectangle function as defined in iBraceweiH 
( 119861) . Then w(R) = a(R) * s(R) which has the form 



w(R) oc e - R2/2c7 ** 



filib/o-al 



dr] 7} e~ v ' /2 I [ 77— I , (19) 

Catm / 
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where Iq is a modified Bessel function of the first kind. We 
have evaluated Equation (T% numerically, and found that as 
long as /?fib/c a tm ;$ 1-5, the convolution of the Gaussian and 
the top-hat will remain approximately Gaussian, where 



CTatm ~ cr atra ^ 1 + % 2 /4 + X 4 /40 

and x = ^fib / Catm- Thus, we will later use 



(20) 



(21) 



for the aperture weighting function in Equation (123b . For 
each lens in our analysis, we take for (T atm the median value 
recorded by the spectroscopic guide cameras during the SDSS 
observations. 

The actual velocity dispersion measured by the observa- 
tions has been effectively luminosity-weighted along the line 
of sight and over the effective spectrometer aperture. This 
averaging can be expressed mathematically as: 



f °° dR R w(R) JZ dZ v{f) (l o].{f) 
J Q c °dRRw(R)J_ZdZ v{r) 



(22) 



where w(R) is the observational aperture weighting function 
defined above in Equation (fJTJ). In Equation ( 1221 the factor 



l-/3yr J takes into account how the radial and tangential 
components of the velocity dispersion tensor project along the 
line of sight. This integral expression for \\^ is analytic 
and, given the above definitions, it becomes 



As 



2 (2^ 



(1+7) 4 Dls 
A(0-/3A(£+2) 



\(a)\(5) 



(C-2/3) 

r(¥) 

3 — S \ 



(23) 



where, again, A(x) is the ratio of gamma functions defined in 
Equation dl5l . Henceforth, we abbreviate ( a 1 n ) as simply 



a\. One can immediately observe from Equation d23l that 
there are degeneracies among a, (3 and 8. The dependence 
upon Oa enters this relation by way of the ratio Ds/D^. 

4. DETERMINATION OF 7 AND f2 A 

As mentioned previously, the general approach to generat- 
ing a constraint is to compare the velocity dispersion from the 
SDSS observations, ctsdss, with the velocity dispersion calcu- 
lated from a galaxy model <r*. By virtue of the analysis in the 
previous section, we have 



:<7*(a,/3,5,6>£;7>^A) 



(24) 



where the semicolon separates the galaxy parameters and ob- 
servables (a, (3, S, 9e) for which we have measured or adopted 
values from the global parameters (7, J7a) which we are seek- 
ing to constrain. 

We will generate the constraint by considering <tsdss an d 
its corresponding uncertainty esdss to be a measurement with 
Gaussian errors such that the probability density for the ob- 
served value of (Tsdss given a "true" value of ct* is 



f(ersDSs|cr*) = 



1 



2-K£ 



■ exp 



SDSS 



(c^-osdss) 



2f 2 

zfc SDSS 



(25) 




Figure 1. Constraint on 7 determined using the method discussed in the text. 
The gray curves represent the posterior PDF for 7 from each lens system. The 
black curve is the joint posterior PDF for all systems, the normalized product 
of the gray curves. (In this plot the joint PDF is scaled by a factor of one-half 
so as to relatively enhance the scale of the individual gray PDF curves.) A 
Gaussian fit to the joint posterior PDF gives 7 = 1.01 ± 0.05. 

Statistical errors on 5, 8e, and the source and lens redshifts 
are negligible in comparison to the velocity dispersion errors 
and the intrinsic variation of the galaxy parameters a and (3, 
and therefore we do not treat them explicitly. 

We determine observational constraints for both 7 and SI a 
in a similar manner and thus the following discussion con- 
cerns our method of constraining a single parameter X £ 
In order to constrain 7 we assume a standard flat 
cosmology with J7a = 0.7. In order to constrain Ha, we as- 
sume that GR is correct and thus 7=1. 

We are interested in calculating the probability of our ob- 
servations having yielded a certain value of galaxy velocity 
dispersion, ctsdss, given a particular assumed value of the pa- 
rameter X. We incorporate the dependence upon S and Oe 
by using the system-specific values determined from the HST 
imaging as discussed in Sections |2] and [3] although in what 
follows we suppress 5 and 9e for notational convenience. Be- 
cause we cannot independently measure a and (3 for individ- 
ual lensing systems, we consider these values to be drawn 
from Gaussian distributions P(a) and P((3) with known mean 
and intrinsic scatter. For a given lens, the probability density 
of interest can be written as 

^sdss|X) = J daP(a) J df3P(f3)P[a SDSS \a4a,(3;X)] . 

(26) 

Assuming a flat prior P(X) for X over a range of interest, a 
universal value for X in all systems, and statistical indepen- 
dence of the measurements of each system, the posterior prob- 
ability for X given the data is then proportional to the product 
over the individual probabilities: 

P(X\ {ctsdss,;}) ocP({<7 S dss,;} \X)P(X) 
53 



:P(X)~[[P(a smS j\X) . (27) 



We illustrate this analysis by selecting a fiducial set of pa- 
rameters. For a and (3 we choose a set of Gaussian distribu- 
tions characterized by: 



(a) = 2.00 ; a a = 0.08 
(/3) = 0.18 ; a (j = 0.13 



(28) 
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The cited 1-er ranges in a and (3 are meant to indicate the 
intrinsic spreads in these quantities rather than uncertainties 
in the mean values. The mean of a is fiducially taken to be 
simply 2.00, the slope of a singular isothermal s phere. This 
is the same value as adopted in the analysis of iGrillo et "all 
(120081) . This is also consistent with the value of (a) = 1.96± 
0.08 determined for th e same sample of SLACS lenses by 
iKoopmans et al.l (120091) via an ensemble aperture-mass anal- 
ysis, which is independent of any dynamical data or model- 
ing and hence does not introduce circularity into our logic 
in proceeding to constraints on 7 and Oa- The scatter in a 
and the distribution for j3 are the same as used in lBolton et al.l 
d2006bl) . and are taken to describe the distribution of mass- 
dynamical properties of the well-studied sample of nearby el- 
liptica l galaxies from lGerhard et all (1200 lb . IKoopmans et al.l 
(120091) present a significantly higher value for the intrinsic 
scatter in a, but argue that this value should be regarded as 
an upper limit given the likelihood of unmodeled systematic 
effects. For comparison with these adopted values of (a) and 
(0), our measured values of S have a mean of (6) = 2.40 and 
a standard deviation as = 0. 1 1 . 

Performing the analysis discussed above for 7, we find the 
resulting posterior probability density shown in Figure \T\ A 
fit to a Gaussian gives 7 = 1.01 ±0.05 (1 a confidence). The 
result is consistent with 7 = 1 and GR. Compared to the pre- 
vious iBolton et al.1 (l2006bl) result, 7 = 0.98 ± 0.07, the sta- 
tistical errors are reduced by a factor of ~ 1 .5 as a result of 
the increase in the number of lensing systems considered (15 
vs. 53). However, as we discuss in the next section, the sys- 
tematic errors may be comparable to, or even dominate over, 
these statistical uncertainties, thereby preventing a large im- 
provement in the precision of 7 simply by increasing the sam- 
ple size of gravitational lens systems. 

We now set 7 = 1 and establish a constraint on Oa- Because 
the calculation, and in particular the ratio Ds/£>ls contained 
in Equation ( 1231 , is not very sensitive (at more than the ~ 10% 
level) to cosmological parameters, we restrict our attention to 
flat universes with Om + Oa = 1. The resulting posterior prob- 
ability density for Oa is shown in Figure [2] The peak occurs 
at Oa = 0.75 and the distribution has a best-fit Gaussian width 
(1 a) of 0.17. Oa = is clearly excluded. This constraint is 
consistent with the results from the cosmic microwave back- 
ground, galaxy clusters, and type la supernovae. The re- 
cent WMAP5 data in combination with type la supernovae 
and baryon acoustic oscillations indicate Oa = 0.726 ± 0.015 
and is also consistent with our assumption of a flat universe, 
giving the constraint -0.0179 < O* < 0.0081 dKomatsu et al.l 
2009). 

One can see from Figure |2]that none of the individual pos- 
terior probability density function (PDF) curves peak: all are 
either strictly increasing or decreasing. This bimodality oc- 
curs because of the relative insensitivity of Ds/I>ls to Oa- 
A quick test shows that a reduction of the errors on ctsdss 
by a factor of ~ 2 would induce a few of the flatter likeli- 
hood curves to peak. However, for most systems, there is 
no value of Oa g [0, 1] such that a*({a),{/3),5;Q\) = ctsdss. 
which leads to the monotonic likelihood curves. 

One may also observe that the posterior PDF curves for two 
of the systems rise or fall more sharply than the others. These 
two lens systems (J0252+0039 and J0737+3216) are among 
the high end of the SLACS sample in terms of lens redshift 
(Zl ~ 0.3) but appear to be otherwise unremarkable. Nei- 
ther lens strongly influences the final result: removing one 
or the other system shifts the peak posterior probability for 
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Figure 2. Constraint on SIa determined as discussed in the text. The gray 
curves represent the individual posterior PDFs of Q\ from each lens system. 
The black curve is the joint posterior PDF, the normalized product of the gray 
curves. A Gaussian fit to the joint PDF gives Qjy = 0.75 ± 0. 17. 
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Figure 3. Sensitivity of the ratio Dls / Os to the cosmological parameter f2 a • 
The gray curves shown are each of the 53 SLACS lenses. The median is 
shown by the thick black curve. All curves are normalized to a value of unity 
at a nominal value of Q/y = 0.7. 

Oa by about three-quarters of a standard deviation, and re- 
moving both leaves the peak of the posterior probability ap- 
proximately unchanged. 

To illustrate the dependence of the ratio D^s/Ds to the cos- 
mological parameter Oa, we show in Figure [3] the ratio as a 
function of the assumed value of Oa for the 53 SLACS lenses. 
The curves are all normalized so that they yield a ratio of unity 
for a nominal value of Oa = 0.7. We see clearly from the 
curves in Figure [3] that the ratio of distances typically varies 
by ±7% and rarely by more than ±10%. The fractional statis- 
tical error in the inferred distance ratios for each system will 
be roughly 10-15%, or twice the fractional error in the mea- 
sured velocity dispersion (Ds/Dis cx a\ ; see Equation (l23l). 
Therefore a determination of Oa from an individual lens sys- 
tem can only be ma de at a "less-than- one-sigma" level. As 
previously noted by IGrillo et all (120081) . many lenses are re- 
quired to obtain a statistically significant determination of Oa- 

5. THE EFFECT OF SYSTEMATICS 

Figure Q] shows the posterior PDF for the post-Newtonian 
parameter 7 using fiducial values for the lens model param- 
eters: (a) = 2.00, a a = 0.08, (0) = 0.18 and a (i = 0.13 (see 
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Figure 4. Derived value of the post-Newtonian parameter, 7, as a function 
of each of the six lens parameters, as defined in the text, centered on their 
fiducial value. As one parameter is varied, the others are held at their fiducial 
values. 

Equation (ESI), the measured slopes of the light profiles, <5, 
and the values of the measured velocity dispersions, (Jsdss- 
The figure provides a clear estimate of the statistical uncer- 
tainty in 7, such that 7 = 1 .01 ± 0.05 (1 a confidence). Simi- 
larly, Figure [2] shows the posterior PDF for the JIa parameter 
using the same fiducial values for the lens model parameters 
as described here for determining 7, which yields a value for 
f^A and its statistical uncertainty, SI a = 0.75 ±0.17. 

However, each of the six lens parameters discussed above 
is susceptible to systematic error in either the adopted mean 
values ((a) and (/?)), the assumed intrinsic widths (<r a and 
073), or in the directly measured parameters (<5 and ctsdss)- It 
is therefore of crucial importance to test the sensitivity of our 
results for 7 and SI a to systematic shifts in the parameters as- 
sociated with the lensing galaxies and to judge whether these 
sensitivities pose immediate and fundamental limitations to 
the precision of our constraints. To this end, we map out the 
dependence of the peak value of 7 and S7a ( e -jL > from results 
that are analogous to those shown in Figures [T] and [2]i on the 
lens parameters. In particular, we vary (a) over the range 
1.90-2.10 and (/3) over the range 0.08-0.28. We explore 
the effect of the assumed intrinsic widths by varying a a over 
the range 0.03-0.18 and ap over the range 0.03-0.23. Next 
we shift all the measured values of dispersion velocity, ctsdss 
by a multiplicative factor, /, taken between 0.95 and 1.05, as 
well as all the measured slopes of the light distribution, S, by 
an additive amount A ranging from -0.10 to 0. 10. As we vary 
each of the these parameters in turn, the other parameters are 
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Figure 5. Derived value of Ha as a function of each of the six lens parame- 
ters, as defined in the text, centered on their fiducial value. As one parameter 
is varied, the others are held at their fiducial values. 



held fixed at their fiducial values. 

The results for the sensitivity of 7 to systematic variations 
in the six lens parameters are shown in Figure[4] In each of the 
six panels the center of the plot is located at 7 = 1 and at the 
fiducial value of the lens parameter being varied: (a) = 2.00, 
(fi) = 0. 18, a a = 0.08, (7^=0.13,/=! .0, and A = 0.0, respec- 
tively. The values of the slopes of these curves (at the fidu- 
cial point) are summarized in Table Q] The analogously de- 
termined sensitivity of SIa to systematic variations in the six 
lens parameters is shown in Figure [5] with the corresponding 
slopes also given in TableQ] The results are closely related to 
those for 7, except that they are typically a factor of ~3 times 
larger in the case of SIa than for 7. This is due to the fact that 
matching the model (cr*\\) (as defined in Equation [23l to the 
measured a is much more sensitively dependent on 7 than it 
is on SI a- Roughly speaking, SIa is about three times more 
difficult to determine precisely than 7 using strong lensing. 

Ingeneral, one can see from Figures H] and [5] and Ta- 
ble QJ that constraints on 7 and SIa are quite sensitive to 
small systematic shifts in the adopted lens-galaxy param- 
eters. For example, a shift in the mean value of the 
slope of the mass profile of 0.04 in (a) leads to a shift 
in 7 of 0.04 x 2.3 ~ 0.09, a change that is nearly twice 
the statistical uncertainty associated with our 7 measure- 
ment. Such a departure of (a) from the value of 2.00, 
which amounts to just one-half of the adopted intrinsic scat- 
ter, is entirely plausible given the variation in this quantity 
as deduced from independent analyses and data sets (e.g., 
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Table 1 

Derivatives with Respect to Lens Parameters 





Y b 


d-y/dY 


dfl A /d¥ 


(a) 


2.00 


2.3 


7.3 


w 


0.18 


0.64 


2.3 




0.08 


0.85 


2.4 




0.13 


-0.03 


0.06 


f 


1.00 


-4.0 


-12 


A 


0.00 


-0.9 


-2.9 



a All lens parameters are dimensionless and are defined in the text. 
b The fiducial value at which the derivatives are evaluated. 



Gerhard et all 120011: IRusin et al.1 12001 IRusin & Kochanekl 
20051: iBolton et al.ll2008bt iKoopmans et al.ll2009l) . More cru- 
cially, this same shift in (a) would produce a change in 
of 0.29, a nearly fatal degradation of our constraint on this 
cosmological parameter. Shifts of comparable magnitude in 
7 and 17a would be caused by a systematic error of 2.5% in 
the measured velocity dispersions, (i.e., / = 1 ± 0.025) which 
is well within the range of uncertainty due to possible mis- 
match between the stellar populations of the lens galaxies and 
the stellar spectrum template set used for the measurement of 

0SDSS- 

Though the effects are a bit less dramatic, similar state- 
ments hold for the sensitivity of the 7 and results to the 
adopted intrinsic mass-slope scatter a a , the mean velocity 
anisotropy (/?}, and the systematic uncertainty in the light- 
profile slope A. With respect to the last quantity, we note 
that although the power-law luminosity profile model that 
we have adopted could be improved upon with a different 
choice such as deVaucouleurs or Sersic, the best-fit models 
from these classes still show significant departures from the 
HST imaging data. Thus the sensitivity of 7 and S^a con- 
straints to systematic errors in the adopted luminosity pro- 
file is unlikely to improve much beyond the level presented 
here. Even massively parametric models for the luminosity 
profile such as the radial B-spline dBolton et al.l2006all2008al) 
would be afflicted by systematic uncertainties related to ellip- 
ticity and deprojection. In addition, other curr ently unmod- 
eled systematic effects such as selection effects dDobler et al.l 
l2008t |Mandelba um et al.ll2009l) and environmental overdensi - 
ties (lAugerl2008l:lTreu et al.l2009tlGuh"riaraes & Sodrel2009l) . 
which appear to contribute only minorly to analyses of the 
structure of SLACS lenses under the assumption of 7 = 1 and 
ACDM, could have significant implications for the type of 
constraints considered in the current work. 

From this analysis, we conclude that the systematics of the 
lens parameters must be controlled to a much tighter degree 
than currently appears possible before a more precise value 
of the post-Newtonian parameter 7 can be determined, and 
before a meaningful independent determination of Oa can be 
derived from strong gravitational lensing. 

6. SUMMARY AND CONCLUSIONS 

In this work we have shown that, for a well-motivated set 
of fiducial lens-galaxy parameters, the current sample of 53 
SLACS gravitational lensing galaxies can constrain either the 
post-Newtonian parameter 7 (on kpc length scales), or the 
cosmological parameter S^a- The statistical errors are such 
that the measurements provide an interesting level of preci- 
sion (~0.05 in 7 and ^0.17 in Ha)- We have paid particular 
attention to quantifying the effects of systematic uncertain- 
ties in the lens parameters on the determination of 7 and Oa- 



These systematic errors are likely to be of equal or greater 
magnitude than the statistical uncertainties, so that the final 
result is largely limited by our imperfect knowledge of the 
galaxy parameters used in the analysis. In particular, the 
derivatives of both 7 and J7a with respect to several impor- 
tant lens parameters (enumerated in §5) have been calculated; 
these are displayed graphically in Figures |4] and [5] and tabu- 
lated in Table Q] 

The larger lens s ample we use as compared to that of 
IBolton et al.l (I2006bl) has enabled us to reduce the statistical 
error on 7 by a factor of ~ 1.5 (from 0.07 to 0.05). However, 
since the uncertainty in 7 due to the systematic errors in the 
lens parameters is likely to be at least as large as the statisti- 
cal errors, simply increasing the lens sample will not improve 
the constraints on 7 until the lens parameters are better under- 
stood. 

Measurements of J7a based on strong lensing involve dif- 
ferent sets of assumptions than other experimental methods, 
and therefore could provide a useful independent estimate of 
this important cosmological parameter. However, our results 
show that in order for strong lensing by galaxies to become 
effective in determining our independent knowledge of 
lens-galaxy mass-dynamical structure will have to improve by 
an order of magnitude. 

Future observations and analysis have some prospect of im- 
proving this situation. High-resolution ground-based spec- 
troscopy with large telescopes and good seeing can reduce 
both the statistical and systematic errors on the measured ve - 
locity dispersions (ICzoske et al.l 120081: iBarnabe et all 120091) . 
and may indeed be able to provide a measurement of the 
average velocity anisotropy for some systems. A more de- 
tailed, spatially resolved dynamical an alysis of nearby galax- 
ies such as from the SA URON survey (lEmsellem et al.ll2004t 
ICappellari et al.l 120061) could be used to derive more accu- 
rate prior distributions on lens-galaxy mass-model parame- 
ters, providing that the differing selection effects between 
SLACS and SAURON samples can be effectively controlled. 
Finally, system-by-system constraints on galaxy mass pro- 
files using higher-order information contained in lensing data 
alone (e.g., lDve & Warrenll2005l: iWillis et aO2006t iDve et al.l 
120081) could be used to remove some of the need for indepen- 
dent priors. 
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